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ABSTRACT 

, Context. Interstellar grain alignment studies are currently experiencing a renaissance due to the development of a new 

, quantitative theory based on Radiative Alignment Torques (RAT). One of the distinguishing predictions of this theory 

is a dependence of the grain alignment efficiency on the relative angle (^) between the magnetic field and the anisotropy 
direction of the radiation field. In an earlier study we found observational evidence for such an effect from observations 
\ of the polarization around the star HD 97300 in the Chamaeleon I cloud. However, due to the large uncertainties in 

O j. the measured visual extinctions, the result was uncertain. 

I ■ Aims. By acquiring explicit spectral classification of the polarization targets, we have sought to perform a more precise 

O ' reanalysis of the existing polarimetry data. 

Methods. We have obtained new spectral types for the stars in our for our polarization sample, which we combine with 
C/) , photometric data from the literature to derive accurate visual extinctions for our sample of background field stars. This 

allows a high accuracy test of the grain alignment efficiency as a function of ^. 

Results. We confirm and improve the measured accuracy of the variability of the grain alignment efficiency with ^, seen 
in the earlier study. We note that the grain temperature (heating) also shows a dependence on ^ which we interpret as 
, a natural effect of the projection of the grain surface to the illuminating radiation source. This dependence also allows 

QQ ■ us to derive an estimate of the fraction of aligned grains in the cloud. 

■ Key words. Radiation mechanisms: general. Techniques: polarimetric, (ISM:) dust, extinction, ISM: magnetic fields 

1. Introduction will likely also provide new probes of the characteristics of 



o 



the dust. 



^ Broadband optical/infra££d in^^tel^^ was ^j^^ ^ standing "text book" explanation for the 

first detected m 1949 (|Hal]! llM9t iHilineiJ 'l949b.a and alignment, via paramagnetic relaxation in rotating grains, 

> was already from the start assumed to be associated with forwar d in the seminal p aper by Davis and 

^ dichroic extinction due to asymmetric dust grains aligned Q^eenstein (Davi s fc GreenstehJUM DG). Over the fol- 

^ . w^kiheiy2ng axis across the m^etic field direction ^^^^^^^^ ^^rious modifications were then proposed 

^ . dmitiMl llMM I Spitzer fc Tukeyl llMa). However, despite ^^j^^^^^ ^1^^ efficiency of this mechanism, including pro 



over 60 years of efforts, the details of the gram align- ppsals for enhancing th e magnetic susceptibility of the ma- 

ment process are still not fully understood. A quantita- ^^^-^j (I.Tones & Spitzerl BoTt IMathisI Il986i) and enhanced 

tive understanding of the physical processes responsible for th e grains due to particle ejections from the 

the gram alignment and hence polarization will allow us ^^^^^^^ dPurcelllllflTflh . Generally, the driving mecha- 

not only a better understanding of interstellar magnetic ■ ^j^^ ^^^j^^^ modifications to 

fields, by e g. putting th e Chandrasekhar-Fermi method ^G mechanism are fixed in the grain's coordinate sys- 

( IChandrasekhar fc Fermi| | l95:^ ) on a more solid footing, but (including the above so called "Purcell rockets" if the 

particle ejection sites are restricted to specific locations on 

the grain surface). 

* Visiting Astronomer, Complejo Astronomico El Leoncito, 

operated under agreement between the Consejo Nacional de Another long-standing possibility to produce gram 

Investigaciones Cientficas y Tecnicas de la Republica Argentina alignment is through mechanical alignment m situations 

and the National Universities of La Plata, Cordoba and San with relative motion between the gas and dust (jGoldlll952l : 

Juan. iLazarian fc Efroimskvl [19961 ) . This flow will tend to cause 
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the grains to align themselves with as small a collisional 
cro ss-section as possible to the flow (iLazarian fc E froimskvl 
Il996). Because most large scale flows in the ISM (including 
winds from hot stars) tend to be ionized, these flows are 
usually constrained to be along the magnetic field lines and 
therefore cause grain alignment with the long axis of the 
grain along the magnetic field direction and polarization 
perpendicular to the magnetic field direction. 

Starting in the 199Q s sev eral authors, including 
ILazarian fc Draind (|l999al lbh and iRobergd (|2QQ4J ) put the 
paramagnetic relaxation paradigm into doubt by showing 
that, due to the internal energy modes of the grains, a dust 
grain will tend to change its orientation in space on rel- 
atively short time scales. Since such a "flip" will cause a 
torque fixed in the grain to drive the grain rotation in the 
opposite direction (in a space-based coordinate system) the 
grain never achieves a significant angular momentum, mak- 
ing paramagnetic relaxation alignment inefficient. 

Parallel to the discovery of these challenges for 
the classical DG model, an alternative theory of in- 
terstellar grain alignment was being worked out by 
Bruce Draine, Alexander Lazarian and their coworkers 



A 



(e.g. iDraine fc Weingartnerll 199 71: ILaz arian fc Hoang''2QQ7') . 
Based on earlv work bv lPolginov fc M itrofanov (1976) they 
showed that an irregular grain with a net helicity will be 
spun up by the differential scattering of the left and right- 
hand circular components in an external light source. Over 
many periods of Lamour precession, the light can then also 
align the grain with their long axis perpendicular to the 
magnetic field, without any contribution from paramag- 
netic relaxation. Because helicity is invariant on reflection, 
this mechanism is not susceptible to the limitations by the 
"thermal flipping" of the grain. The sole environmental 
requirement for alignment by this mechanism is that the 
radiation field be anisotropic. In the interstellar medium, 
this is almost always the case. Hence, while the physical 
mechanisms responsible for the grain alignment are quite 
distinct, to first order extended DG alignment and radia- 
tive alignment provide the same observational prediction; 
namely polarization parallel to the projected magnetic field 
direction. 

Radiative Alignment Torque (RAT) theory has matured 
over the last decade and is now prov iding a number o f spe- 
cific, testable, predictions (Lazaria n fc HoangI |2QQ7[ ). One 
such prediction, which is especially attractive for probing 
the vahdity of the theory, since it is unique to RAT theory, is 
that the alignment efficiency, but not the polarization angle, 
should vary with the angle between the magnetic field and 
the radiation field anisotropy (^; see Lazarian fc HoangI 
|2QQ7[ ). One way to observationally probe this prediction is 
to find dust grains for which the dominant radiation source 
is a localized source (a star) and measure the alignment 
efficiency for background stars projected at different posi- 
tion angles around the source star. Figure 1 illustrates the 
geometry of the situation. 

We found such a combination of star and nearby dust in 
t he C hamaeleon I cloud and showed in Andersson fc Potte^ 
(j2010j) , combin ing new polarim etry with the multi-band po- 
larimetry from IWhittetl (|1992 ^. that the grain alignment is 
indeed enhanced at small projected distances to the central 
source star (in this case HD 97300). A preliminary study of 
th e angle dependence o f the a lignment was also carried out 
in lAndersson fc Potterl (|2010f ) , but was limited by the fact 
that, for most of our field stars, the spectral classification 
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Fig. 1. The geometry used to probe the effects of the rel- 
ative orientation of the magnetic field and radiation field 
anisotropy is illustrated. In the two panels, the dashed lines 
illustrate an idealized magnetic field, (a; side view) A bright 
source is located close to and in front of an underlying in- 
terstellar cloud surface, such that its radiation locally dom- 
inates over the diffuse radiation field, (b; en face view) By 
probing the polarization of background field stars we can 
probe the effects of angle dependent Radiative Alignment 
Torque (RAT) gr ain alignment (Adapted, w ith permission 
of the A AS, from lAndersson fc Potterll2010D 



and therefore the determinations of the visual extinctions 
on each line of sight, had to be accomplished using near in- 
frared (NIR; 2MASS) photometry alone. Because the small 
color excesses in the NIR and uncertainties in the reddening 
vector slope, the spectral classifications typically had uncer- 
tainties of a full spectral class, translating to uncertainties 
on the visual extinctions of ~1 mag. Despite these limita- 
tions, we detected a relative grain alignment enhancement 
for position angles where the magnetic field and radiation 
field anisotropy were parallel, in agreement with the pre- 
dictions of RAT theory. 

In order to test whether the large uncertainties 
in the derived visual extinctions confused the analy- 
sis, we h ere present a reanalysis of the polarization 
data from 'Andersson fc Potterl (|2010l ), using new explicit 
spectral classifications of the field star acquired at the 
the "Complejo Astronomico de Leoncito" (CASLEO) in 
Argentina. 
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2. Observations, data reduction and analysis 

We u sed the REOSC spectrometer (|PintadQ fc AdelmanI 



11996^ on the 2.15 m telescope at the CASLEO on the 
nights of 2011, March 4-6. The spectrometer was used with 
the 300 lines/mm grating, producing a measure wavelength 
coverage of AA3631-7128 A with a spectral resolution of 
AA=3.4 A/pixel. The detector used was a thinned TEK 
1024x1024 CCD with 24 /im pixels. Wavelength calibra- 
tion was achieved using exposures of Thorium- Argon hol- 
low cathode lamps. No attempt was made to achieve pho- 
tometric calibration of the data as the sky conditions did 
not warrant it. 

The data reduction was achieved using standard pro- 
cedures and routines in the IRAF environment. After 
bias and flat-field corrections, the two-dimensional spec- 
tra were traced and integrated across the dispersion di- 
rection and extracted into one-dimensional form. The ex- 
tracted spectra were then wavelength calibrated, normal- 
ized and compared , intera ctively, to the standard sequence 
from iJacobv et al.l (|1984| ) . The classifications were redone 
independently several times, by two of the authors. The 
resulting spectral classifications are listed in Table 1. Also 
listed i n Table 1 are the derived colour excesses, based on 
Tvch o (lH0g et al.l l2000l ) photometry and intrinsic colours 
fromlCoxl (|2000[ ). In a small number of cases Tycho pho- 
tometry is not available, a nd in those cases we have in- 
stead used NOMAD (Zac harias et al.l [2005) photometry. 
Visual extinctions (Ay) were calculated using total-to- 
selective extinction {Rv) values derived, from Ay = Ry x 
Eb-v, using the relation: R y = 1.1 x Ey-x / Eb-v from 
IWhittet & van Breda (1978)^ 

As discussed in A ndersson fc Potter! (|2007f ), the wave- 
length of the maximum of the polarization curve (Amax) 
is a sensitive probe of the grain alignment and one that 
is immune to the magnetic field topology along the line of 
sight. The absolute value of Amax depends on the over- all 
grain size distributions and likely on the color and intensity 
of the local diffuse interstellar rad i ation field (ISRF), but 
as shown by 'Andersso n fc Potterl (|2007l ) the slope in the 
Amax vs. Ay relation is universal. Hence we can measure 
the relative grain alignment enhancement in a region by 
deriving the average relationship between Amax and Ay for 
field stars behind the cloud in question and then probing 
for localized deviations from this average relationship. We 
did this in Andersson fc Potter (2007, 2010) and found that 
in Chamaeleon I, but beyond the region around HD 97300 
(where its radiation field dominates that of the diffuse ra- 
diation), a linear relationship is found: 

(Amax(Ay)) = (0.527 ± 0.006) + (0.033 ± 0.003) x Ay. (1) 



3. Results and discussion 

3.1. Stellar classification 

The uncertainties for the spectral classifications are based 
on estimates from the two independent classifiers. Usually 
the results by the individual classifiers agreed to within the 
assigned uncertainties. The luminosity class assignments 
are less certain. For early-type star (up to F) we assume 
that the stars are on the Main Sequence. For later stellar 
types (G and beyond), we used the k nown distance to th e 
Chamaeleon I cloud [d = 150 ± 15 pc. lWhittet et al.lll997[ ). 
spectroscopic parallaxes for luminosity classes III and V, 



and the measured colour excesses, derived for each lumi- 
nosity class, to set the luminosity class. The restriction to 
these two luminosity classes could be a source of some er- 
ror, in a few cases where luminosity classes IV (or I) might 
have been more appropriate. Specifically, TYC 9410-1931-1 
and TYC 9414-0046-1 have been assigned a spectral class 
of G4 III. However, this spectral/luminosity class is located 
in the Herzsprung gap of the HR diagram where the num- 
ber of stars is quite small. As noted below some stars show 
colour excess ratios that would indicate that they are part 
of binaries or show the effects of peculiar circumstellar ma- 
terial, or other effects yielding non-standard colours. The 
former of these two stars falls in this category. 

For five stars in our sample the colour excess Eb-v is 
consistent with zero and we have excluded these stars from 
the subsequent analysis. Similarly, for six stars the colour 
excess ratio Ey-K / Es-y is significantly smaller than the 
nominal value 2.74, assuming a st andard in terstellar extinc- 
tion curve and a Ry value of 3.1 (|Coxll200Q I). After checking 
the spectral classifications for these stars we conclude that 
the discrepancy is in the measured colours and may indi- 
cate binarity, or other non-standard spectral behaviour. We 
have therefore excluded also these stars from the analysis 
of the grain alignment variations. 

While a systematic increase in Ry is expected for grain 
growth, no clear correlation is seen between Ry and column 
density (e.g. /loo)- This is likely because the total column 
densities in our sample are still fairly small. 

3.2. Grain alignment 

Figure [2] shows the relative grain alignment efficiencies for 
the stars in our sample as measured by the above technique 
in filled (black) diamonds as a function of ^ (see Figure [T]): 



^^max ^max 



(A„,ax(^V)>. 



(2) 



The FIR ratio leo/Iioo is also shown, in open (red) dia- 
monds. The FIR dat a here are from the IRIS reproce ssing 
of the IRA S data (.Miville-Deschenes fc Lagachell2005l ). For 
details see I Andersson fc Potterl (|2007[ ). 

In Figure [2] we have over-plotted the best fit of the func- 
tion AAmax = A -\- B X cos(^ — ^q) whcrc ^ is the posi- 
tion angle from the star, relative to the average magnetic 
field direction in the area, as de termi n ed from the polariza- 
tion maps by iMcGregor et al.l (|1994| ) ; lAndersson fc Potterl 
(2010) and is the symmetry angle of the function (in 
this case the minimum). We have chosen a simple cosine 
relation here, since the exact theoretical functional form 
for the grain alignment depends on several unknown pa- 
rameters of the grains, the radiation field and the structure 
of the cloud along the line of sight. With the new and sig- 
nificantly improved visual extinctions for our background 
stars, we find a statistically significant depression in AAmax 
^0 close t o 0, in agreement with theoretical predictions. 
An F-test (jLup ton 1993) yields a greater than 99% prob- 
ability that the two additional parameters in the cosine 
function are statistically warranted (as compare to a sim- 
ple average value). The parameters in the current best fit 
are all within the mutual uncertainties of the best fit in 
lAndersson fc Potterl (|2010t ). We now find a 9.5cr deviation 
from a null result for the amplitude of the grain alignment 
enhancement. 
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Fig. 2. The variations in grain alignment efficiency, as 
probed by the location of the peak in the polarization curve, 
is shown as a function of the relative angle between the mag- 
netic field and radiation field anisotropy: ^ (filled, black, 
diamonds. Also shown are the variations in the FIR ratio 
^6o/^ioo (open, red, diamonds). Both show a variation with 
^. The change in the grain temperature can be understood 
in terms of the surface area facing the illuminating star, 
while the variations in Amax are consistent with predictions 
by RAT theory. 

The fractional polarization (pmax/^y; Fig [3]) for our 
combined sample shows a power-law slope b = —0.39 ±0.09 
consistent, within 1.2a mutual uncertainty, from the results 
in^hittet et al.l ([2008): b = -0.52 ± 0.07. 

As discussed in Anderss on fc Potterl (j2010|) , while we do 
detect an alignment efficiency enhancement in the direction 
of the magnetic field, relative to the star, we do not see any 
evidence for a rotation of the polarization angles in this 
direction, as would b e exp ected for mechanical alignment 
(see Figure 9 of I Andersson fc Potterl l2010l ). 




Av [mag] 

Fig. 3. The fractional polarization (pmax/^y) for the sam- 
ple is plotted as a function of the visual extinction. The 
power-law exponent is consistent, within the mutual error- 
bars, with the one found for deeper lines of sight by 
IWhittet et J](|2008D . 
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3.3. Grain heating 

Figure [2] also shows a systematic variation in the /eo/Aoo 
ratio as a function of ^. To model this we envision the 
alignable grains as a system of simple parallelepipeds with 
two large sides of length a, and a smaller side of length b 
(like a pizza box; Figure H]). The radiative heating of these 
grains around the star then takes the form: 

/ oc ^a^[cos\^) + -sin2(^)], (3) 

where Iq is the effective intensity of the central star and 
r is the distance between the star an dust grain. Figure |4] 
illustrates the model geometry. We also need to allow for 
symmetrical grains that will not show a d ependence on the 
angle ^. As shown bv .Desert et al.l (|l990l) the /eo/^ioo ra- 
tio shows a roughly linear response to the radiation field 
strength in the range 1-5 times the local interstellar radia- 
tion field. Hence, as a simple first order, heuristical, model 
we can fit the FIR colour temperature to the function: 

^(M/) = ci[(l - e) + e(cos2(M/) + (b/a)sin^m], (4) 

^100 



Fig. 4. The heating of asymmetrical grains, aligned with 
the magnetic field (thin, gray, lines) will vary with the an- 
gle ^. If only a fraction of the grains are alignable, we can 
use the variations in the dust temperature to estimate the 
fraction of aligned grains or the axes ratio of these asym- 
metric grains 



where e is the fraction of aligned, asymmetrical grains and 
b/a is the axes ratio of the asymmetrical grains. 

The two parameters e and (b/a) are degenetrate in the 
fits. We therefore ran a number of fits using e as a free pa- 
rameter with {b/a)~^ set to y^2, 2, 4 and 6. For these choices 
of {b/a) we find e=(0.58±0.04), (0.34±0.02), (0.23±0.02) 
and (0.20±0.02). Figure [2] shows the best fit for (6/a)-^=4. 

iKim fc Martini (|1995) present, in their Table 1, the 
polarization-to-extinction ratio for perfectly aligned oblate 
spheroids with the axes ratios we used in our fits. They 
also quote the maximum of the measured value of this ra- 
tio as Pmax/T=0.028. If we use the ratio of the observed to 
calculated value in their paper as signifying the required 
dilution of the aligned oblate spheroidals (by unaligned. 
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or fully symmetrical grains), we find that the fractions of 
aligned grains they predict are 0.56, 0.31, 0.18 and 0.17, 
very close to the best fit values in our model. 

The fitting constant ci in equation S] contains the dis- 
tance from the illuminating star to each line of sight. Our 
sample is not large enough to conclusively show this depen- 
dence for different distance annuli, but splitting the sample 
in two and fitting the targets nearest to HD 97300 does 
produce an increased value of ci , compared both to the full 
sample and the more distant half-sample. 

Repeating this analysis on the full maps of IRAS data 
in the Chamaeleon region produces consistent results al- 
beit with larger scatter in the FIR ratio as a function of 
position angle (in particular, for the annuli d=0.3-0A deg 
and d=0A-0.5 deg, the best fit values for e agree, within 
the error bars, with the ones found for our polar imetry star 
sample). Further analysis, in particular repeating the ex- 
periment for other cloud regions dominated by a near-by 
star, will be needed to ascertain whether this effect is real 
or whether it is a statistical aberration in the current data 
set. 



4. Conclusions 

With the improved o bservational data, we confirm the pre- 
liminary result from lAndersson fc Potterl (|2010f ) that the 
grain alignment seems to depend on ^, as predicted by RAT 
theory. This provides an important observational constraint 
on, and in this case support for, the currently best devel- 
oped theory for interstellar grain alignment. With further 
supportive observational tests of the theory, the second- 
longest standing myster^Q of interstellar medium astro- 
physics may be within reach of resolution. 

We also find a dependence on ^ in the leo/hoo ra- 
tio. Since the aligned grains causing the polarization have 
their major axis perpendicular to the magnetic field, this 
enhancement in the FIR ratio should be expected simply 
as a projection effect of the grain surface to the radiation 
from the central star. Our best fit models of the FIR ratio 
find that the fraction of aligned grains is very close to that 
found for theo retical models o f the extinction and polariza- 
tion curves bv lKim fc Martini (p^995). 

It is possible that an underlying corrugation of the cloud 
surface could partly be responsible for the ^ dependence. 
This could happen if a ridge centered on HD 97300 is ori- 
ented parallel to the average magnetic field. However, the 
low estimated magnetic field of the Chamaeleon I cloud, 
B\\ = 3 ± 4 /iG (|Bourke et all 120011 ) makes this unlikely. 
Repeating the experiment in other regions where a local- 
ized radiation source dominates the grain illumination will 
allow such caveats to be addressed. 



^ The longest-standing mystery of ISM astrophysics is the na- 
ture of_th^_caxriersof_^ Diffuse Interstellar Bands, detected 
fsee iFriedman et al.ll20Tll ) in 1919 ( Heger 1922) and identified 
as being of interstellar origin in 1936 (f Merrill. 193d ) 
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Table 1. Spectral types and photometry for stars in the region. 
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The number in parenthesis gives the estimated uncertainty of the spectral class, in sub-classes; ^ Source of spectral 
classification: A - This work, B - Michigan Spectral Survey (|Houk Swift 1999), C - SIMBAD database; ''^^ In a few cases 
unphysical results were obtained for Rv- These likely indicate non-standard spectral classes. Those values have been omitted and 
marked by — and the resulting visual extinctions were not used in the analysis. 
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